INVERSE SPECTRAL PROBLEMS FOR STURM-LIOUVILLE 
OPERATORS WITH SINGULAR POTENTIALS, II. 
RECONSTRUCTION BY TWO SPECTRA* 

ROSTYSLAV O. HRYNIV AND YAROSLAV V. MYKYTYUK 

Abstract. We solve the inverse spectral problem of recovering the singular poten- 
tials q G W 2 -1 (0,1) of St urm-Liouville operators by two spectra. The reconstruction 
algorithm is presented and necessary and sufficient conditions on two sequences to be 
spectral data for Sturm-Liouville operators under consideration are given. 



I. Introduction 

Suppose that q is a real- valued distribution from the space W% (0, 1). A Sturm- 
Liouville operator T in a Hilbert space 7i := L 2 (0, 1) corresponding formally to the 
differential expression 



can be defined as follows |18j . We take a distributional primitive a G 7i of q, put 

l a {u) := —{v! — cru)' — ou 

for absolutely continuous functions u whose quasi-derivative := u' — au is again 
absolutely continuous, and then set 

Tu := l a (u) 

for suitable u (see Section |21 for rigorous definitions). When considered on this natural 
domain subject to the boundary conditions 

(1.2) u [1] (0) - Hu(0) = 0, m [1] (1) +hu(l) = 0, H,h E MU {oo}, 

the operator T = T(cr, H, h) becomes [TTHIH] selfadjoint, bounded below, and possesses 
a simple discrete spectrum accumulating at +oo. (If H — oo or h — oo, then the 
corresponding boundary condition is understood as a Dirichlet one.) 

We observe that among the singular potentials included in this scheme there are, e.g., 
the Dirac 5-potentials and Coulomb l/a;-like potentials that have been widely used in 
quantum mechanics to model particle interactions of various types (see, e.g., PQI2])- It 
is also well known that singular potentials of this kind usually do not produce a single 
Sturm-Liouville operator. Still for many reasons T(cr, H, h) can be regarded as a natural 
operator associated with differential expression ()1.1|) and boundary conditions (jl.2|) . 
For instance, l a {u) = —u" + qu in the sense of distributions and hence for regular (i.e., 
locally integrable) potentials q the above definition of T(a, H, h) coincides with the 
classical one. Also T(cx, H, h) depends continuously in the uniform resolvent sense on 
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cr G TC, which allows us to regard T(a, H, h) for singular q = a' as a limit of regular 
Sturm-Liouville operators. 

Sturm-Liouville and Schrodinger operators with distributional potentials from the 
space W2 -1 (0, 1) and W^Vum/C^) respectively have been shown to possess many prop- 
erties similar to those for operators with regular potentials, see, e. g., |§1 ITTH IT7| ITS]. 
In particular, just as in the regular case, the potential q = a' and boundary condi- 
tions (jl.2j) of the Sturm-Liouville operator T(a, H, h) can be recovered via the cor- 
responding spectral data, the sequences of eigenvalues and so-called norming con- 
stants [TTj . 

The main aim of the present paper is to treat the following inverse spectral prob- 
lem. Suppose that {A^} is the spectrum of the operator T(cr, H, hi) and {^1} is the 
spectrum of the operator T(cr, H, h 2 ) with h 2 ^ hi. Is it possible to recover a, H, 
hi, and h 2 from the two given spectra? More generally, the task is to find necessary 
and sufficient conditions for two sequences (A^) and (fj^) so that they could serve as 
spectra of Sturm-Liouville operators on (0, 1) with some potential q = a' G W 2 ~ 1 (0, 1) 
and boundary conditions (jl.2|) for two different values of h and, then, to present an 
algorithm recovering these two operators (i.e., an algorithm determining a and the 
boundary conditions). 

For the case of a regular (i.e., locally integrable) potential q the above problem was 
treated by Levitan and Gasymov [UJ US] when hi and h 2 are finite and by Mar- 
CHENKO [To^ Ch. 3.4] when one of hi, h 2 is infinite. Earlier, Borg jl] proved that two 
such spectra determine the regular potential uniquely. 

Observe that T(cr + h',H — h', h + h') = T(cr, H, h) for any h' G C, so that it 
is impossible to recover a,H,hi, and h 2 uniquely; however, we shall show that the 
potential q = o' and the very operators T(a, H, hi) and T(cx, H, h 2 ) are determined 
uniquely by the two spectra. The reconstruction procedure consists in reduction of 
the problem to recovering the potential q = a' and the boundary conditions based on 
the spectrum and the sequence of so-called norming constants. The latter problem 
has been completely solved in our paper [TTj, and this allows us to give a complete 
description of the set of spectral data and to develop a reconstruction algorithm for 
the inverse spectral problem under consideration. 

The organization of the paper is as follows. In the next two sections we restrict 
ourselves to the case where H = oo and hi,h 2 G M. In Section |2] refined eigenvalue 
asymptotics is found and some other necessary conditions on spectral data are es- 
tablished. These results are then used in Section 01 to determine the set of norming 
constants and completely solve the inverse spectral problem. The case where one of hi 
and h 2 is infinite is treated in Section^ and in Section |S] we comment on the changes 
to be made if H is finite. Finally, Appendix El contains some facts about Riesz bases 
of sines and cosines in L 2 (0, 1) that are frequently used throughout the paper. 



2. Spectral asymptotics 

In this section (and until SectionlH) we shall consider the case of the Dirichlet bound- 
ary condition at the point x = and the boundary conditions of the third type at the 
point x = 1 (i. e., to the case where H = oo and hi,h 2 G M). We start with the precise 
definition of the Sturm-Liouville operators under consideration. 

Suppose that q is a real- valued distribution from the class W 2 l {Q, 1) and a G H is any 
of its real-valued distributional primitives. For h G K we denote by T a ^ = T(a, oo, h) 
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the operator in 7i given by 

T a>h u = l a (u) := -(u [1] )' - au' 

and the boundary conditions 

u(0) = 0, u [1] (1) + hu{l) = 0, heR 

(we recall that vP^ stands for the quasi-derivative u' — au of it). More precisely, the 
domain of T a>h equals 

D(T^) = {ue Wl[0, 1] I u [1] G Wl[0, 1], Z a («) G H, w(0) = 0, + /m(l) = 0}. 

It follows from ^H] that so defined is a selfadjoint operator with simple discrete 
spectrum accumulating at +00. The results of [TH] also imply that the eigenvalues A 2 
of T ay h, when ordered increasingly, obey the asymptotic relation A n = n(n — 1/2) + A n 
with an £ 2 -sequence (A n ). However, we shall need a more precise form of A n , and shall 
derive it next. 

To start with, we introduce an operator T a defined by T G u = l a (u) on the domain 
X3(T CT ) = {ue Wl[0, 1] I u [1] G Wl[Q, 1], l a {u) G H, f(0) = 0}. 

In other words, T a is a one-dimensional extension of T a ^ discarding the boundary 
condition at the point x — 1. We invoke now the following results of on the 
operator T a . Just as in the classical situation with regular potential q the operators T± a 
turn out to be similar to the potential- free operator To; the similarity is performed by 
the transformation operators I + K^, where are integral Volterra operators of 
Hilbert-Schmidt class given by 



Kfu{x) = [ k±(x,t)u(t)dt. 
Jo 



It follows that any function u from D(T a ) has the form u = (I+K^)v for some function 
v G 2D(T ) (i.e., for some v G Wf [0, 1] satisfying the boundary condition v(0) = 0) and 
that T a (I + K£)v = —(I + K+)v". Moreover, there exists a Hilbert-Schmidt kernel 
r ff (x, t) such that 

[(I + Kt)vf\x) = [(/ + K-)v']{x) + / r a (x,t)v(t) dt 







for every v G S3 (To). Also, for every fixed x G [0, 1] the functions k„(x, ■) and r a (x, 
belong to 7i. 

Due to the above similarity of T a and To, for any nonzero A G C the function 



u(x, A) := sin Aa; + / k£(x, t) sin At dt 
Jo 



is an eigenfunction of the operator T a corresponding to the eigenvalue A 2 G C. If in 
addition u(x, A) satisfies the boundary condition w^(l) + hu(l) = 0, then u(x, A) is also 
an eigenfunction of the operator T a ^ corresponding to the eigenvalue A 2 . Therefore the 
nonzero spectrum of T a h consists of the squared nonzero A-solutions of the equation 



(2.1) 



f 1 1 f 1 

cosA-l- / k a (1, t) cos At dt + — / r CT (l, t) sin At dt 

Jo A J 



h ■ ^ 

+ — sin A + 
A 



h r 1 

- / k+{l, t) sin Xtdt = 0. 
A Jo 



4 



R. O. HRYNIV, YA. V. MYKYTYUK 



Recall that the operator T a ^ is bounded below and hence after addition to q a suitable 
constant T a ^ becomes positive. Therefore without loss of generality we may assume 
that all zeros of equation (J2.1|) are real. Observe also that due to the symmetry we 
may consider only the positive zeros. The asymptotics of these zeros is given in the 
following theorem. 

Theorem 2.1. Suppose that a e H and h G K. are such that the operator is 
positive and denote by X\ < X\ < ■ ■ ■ the eigenvalues of T a ^. Then the numbers X n 
satisfy the relation 

X n = n(n - 1/2) + — + 5 n (h) + 7n , 
Tin 

in which the sequences (5 n (h))™ and ( n / n )^=x belong to l\ and £2 respectively with 7 n 
independent of h. 

Proof. Recall that X n — > 00 as n — > 00 and that the functions k^(l,t) and r a (l,t) 
belong to Ti. Therefore (|2.1j) and the Riemann lemma imply that 

cosA n = — / k~(l,t) cos\ n tdt + 0(1/X n ) = o(l) 
Jo 

as n —* 00. By the standard Rouche-type arguments (see details, e.g., in [TBI Ch. 1.3]) 
we conclude that A n = 7i{n — 1/2) + A n with A n — > as n — > 00. As a result (see 
AppendixEJ), the system {sin X n x} turns out to be a Riesz basis of TC and the sequences 
(CO*! and «)~ 1 with 

c n := / k„(l, t) sin X n t dt, c' n := / r CT (l, t) sin A„t dt 
Jo Jo 

belong to £2- Relation ()2.1|) now implies that 

Sin X n = (-l)^ 1 / k~(l, t) COS X n t dt + + (_l)n+l fec n + C n = 

Jo A n A ra 

= Tn + — + <^n + C 

7rn 



where 



»i 

7n :=(-l) n+1 / k-(l,t)coan(n-l/2)tdt, 
Jo 

(2.2) 5' n ,:= (-l) n+1 / feJ(l,t)[cosAn*-cos7r(n-l/2)t]dt, 

frcosA„. /i n „+i^c n + c n 

5 n := x — + I" 1 ) 



A n vrn X n 

Next we observe that the systems {cos 71(72 — l/2)t} and {cosA„t} form Riesz bases 
of H, so that (7„) and (8' n ) belong to £ 2 - It follows that (sinA„) G £ 2 and therefore 
(A n ) G £2- Lemma lA.21 now implies that the numbers S' n form an ^-sequence, and 
simple arguments justify the same statement about 5". 

Thus we have shown that 



. ~ h . 

sin X n = 7„ H h d n + 6: 

Tin 
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with the sequences (7„) from £ 2 and (S' n ) and (5'^) from l\. Applying arcsin to both 
parts of this equality, we arrive at 

~ ft . 

An = In H h n 

TTTl 

with some £i-sequence (S n ), and the theorem is proved. □ 

It was shown in ^Tj that any sequence (X n ) of pairwise distinct positive numbers 
satisfying the relation A„ = ir(n — 1/2) + A„ with an £ 2 -sequence (A„) is the spectrum of 
some Sturm-Liouville operator T a> h with potential q = a' from Wjf 1 (0, 1), the Dirichlet 
boundary condition at x — 0, and third type boundary condition (jl.2j) at x — 1 with 
some /igR. Suppose that (fi n ) is the spectrum of a Sturm-Liouville operator with the 
same potential q, the Dirichlet boundary condition at x — 0, and third type boundary 
condition (jl.2|) at x = 1 with a different ft, e R. We ask whether there exists any 
relation between the two spectra. 

The first restriction, just as in the regular case, is that the spectra should interlace. 

Lemma 2.2. Suppose that a &H is real valued and that sequences {X n } and {fi n } are 
the spectra of Sturm-Liouville operators T cr> h 1 and T a ^ 2 with distinct hi, h 2 G M. Then 
these sequences interlace. 

Proof. Denote by u = u(x, A) a solution of equation l a (u) = \u satisfying the boundary 
condition u(0) = 0. We recall that u being a solution of l a (u) = Xu means that 

and hence u enjoys the standard uniqueness properties of solutions to first order differ- 
ential systems with regular (i.e. locally integrable) coefficients. The numbers A n and \i n 
are then solutions of the equations A)+ft,iu(l, A) = and u,M(l, X) + h 2 u(l, A) = 

respectively. 

We introduce a continuous function 9(x, A) through cot6*(x, A) = w^/u. After dif- 
ferentiating both sides of this identity in x and using (|2.3J) we get 

9' —au'u — u^u' , . „ , 9 

— ttt; = 5 A = "( cot 6 + °) ~ A > 

sin 9 u z 

or 

9' = X sin 2 9 + (cos 9 + o- sin 9) 2 . 

It follows from [3J Ch. 8.4] or proof of Theorem XI. 3.1] that the function 9(1, X) 
is strictly increasing in A; hence solutions of the equations cot 9(1, A) = —hi and 
cot 9(1, A) = — ft 2 interlace, and the proof is complete. □ 

The next restriction concerns asymptotics of A n and [i n . 

Lemma 2.3. Suppose that (X n ) and (fi n ) are spectra of Sturm-Liouville operators T^ hl 
and T a ^ 2 with real-valued a ETC and hi, h 2 £ R- Then there exists an £ 2 -sequence (z/„) 
such that 

(2.4) A n - fi n = 1 • 

Tin n 

Proof. Put A„ := A„ — ix(n — 1/2) and /!„ := \i n — n(n — 1/2) so that A„, — /i„ = A„ — /!„. 
Recall that 

sin A„ - sin/2 n = hzJ^l + 5 > ( hl ) _ 5 > (h 2 ) + ^(^) _ S '^h 2 ) 
nn 
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with the quantities 5' n (h),5"(h) of (|2.2j) . Relation 

$'n(hi) ~ S' n (h2) = (-l) n+1 / fe;(l,t)[cosAj-cos//^]^ 

Jo 

and Lemma lA.21 show that 6' n (hi) — b~' n {h 2 ) = (A n — /i n )z/ for some £ 2 -sequence (i/). 
Also there exists an £ 2 -sequence (i^) such that — ^"(^2) = z/ n/ n - Finally, using 

the relations 

sinA n - sin/2„ = (A n - /!„) cosz/ n , 

where i? n are points between A n and Jl n (so that (1 — cosz/ n ) is an £ 2 -sequence), and 
combining the above relations we arrive at 

/T ~ \ / ~ /\ hi — h 2 v'n 

(A n - /i n )(cos^ n - z/J = 1 , 

im n 

which implies (12. 4j) . □ 
Corollary 2.4. Under the above assumptions, lim„_ >00 (A^ — //^) = 2(/ii — /i 2 ). 



3. Reduction to the inverse spectral problem by one spectrum and 

norming constants 

Suppose that a G H is real-valued and that {A 2 } and {/i 2 } are eigenvalues of the 
operators T a>hl and T a>fl2 respectively introduced in the previous section. We shall show 
how the problem of recovering er, hi, and h 2 based on these spectra can be reduced 
to the problem of recovering a and hi based on the spectrum {A 2 } of T a hl and the 
so-called norming constants {a n } defined below. This latter problem for the class of 
Sturm-Liouville operators with singular potentials from W 2 ~ 1 (0, 1) is completely solved 
in [TT]; see also [TSl EHJ for the regular case of integrable potentials. 

Denote by Ui(x, A) and u 2 (x, A) solutions of the equation l a {u) = X 2 u satisfying the 
initial conditions Ui(l, A) = u 2 (l, A) = 1 and Uj(l, A) = —hj, j = 1, 2. Then U\(x, A n ) 
are the eigenfunctions of the operator T cr / ll corresponding to the eigenvalues A 2 , and 
we put 

a n := 2 / \ui(x, X n )\ 2 dx. 



Our next aim is to show that a n can be expressed in terms of the spectra {A 2 , only. 

Set $j(A) := Uj(0,X), j = 1,2; then zeros of $1 and $ 2 are precisely the numbers 
±A n and ±/i n respectively. We show that $1 and $ 2 uniquely determine a n and are 
uniquely determined by their zeros. 

Lemma 3.1. The norming constants a n satisfy the following equality: 

hi-h 2 $i(A n ) 



(3.1) a, 



An $ 2 (A n ) 



Proof. The proof is rather standard and the details can be found, e.g., in [T3], so we 
shall only sketch its main points here. 

Put m(A) := — $ 2 (A)/$i(A); then the function f(x,X) := u 2 (x, A) + m(X)ui(x, A) 
for all A G C not in the spectrum of T CTj / ll satisfies the equation l a (f) = A 2 / and the 
boundary condition /(0, A) = 0. In other words, f(x,X) is an eigenfunction of the 



RECONSTRUCTION BY TWO SPECTRA 



operator T a corresponding to the eigenvalue A 2 . Since Ui(x, X n ) is an eigenfunction of 
the operator T a corresponding to the eigenvalue A 2 , we have 

(A 2 ~ A 2 ) / A)mi(x, Xn) dx = [T a f, Mi) - (/, T a Ui) 



= A)m(l, A) + f(l, A)«W(1, A n ) = h 2 -h 1 . 

On the other hand, 



(A 2 - A 2 ) f 
Jo 



-1 

f(x, X)ui(x, X n ) dx = (A 2 - A 2 ) / u 2 (x,X)u l (x,X n ) dx 



,2 x2^(A) ^ 



(A - A n ) ^ . J ui(x, X)ui(x, X n )dx, 



and after combining the two relations and letting A — > A n we arrive at 

_ h l -h 2 $;(A ra ) 

A n $ 2 (A n )' 

The proof is complete. □ 

The following statement has appeared in many variants in numerous sources, but we 
include its proof here for the sake of completeness. 

Lemma 3.2. In order that a function f(X) admit the representation 

(3.2) /(A)=cosA+/ g(t) cos Xtdt 

Jo 

with an TC-function g, it is necessary and sufficient that 

(33) fW-U^Tjw- 

where fj. := ir(k — 1/2) + fk and is an t 2 -sequence. 

Proof. Necessity. The function / of ()3.2|) is an even entire function of order 1, and the 
standard Rouche-type arguments (see the analysis of Section |2J show that zeros ±fk 
of / have the asymptotics fk = Tt{k — 1/2) + f k with (f k ) e On the other hand, up 
to a scalar factor C, the function / is recovered from its zeros as 

00 \ 2 00 / x2 \2\ 

where 

" 7r 2 (yfc - 1/2) 2 



2 

/,• _ 1 J fc 



ft 

To determine Ci, we observe that, in view of (|3.2|h lim^oo f(iv)/ cosiu = 1 and that 

x ^ (vr 2 (fc-l/2) 2 -A 2 ) 

(3 - 5) cosA = n ^,^2 ; 

comparing now (|3.4|) and (|3.5|) . we conclude that Ci = 1 and necessity is justified. 

Sufficiency. Suppose now that / has a representation of the form (|3.3|) . in which 
fk = 7r(k — 1/2) + fk with (fk) G £2- We assume first that the zeros fk are pairwise 
distinct. Then due to the asymptotics of fk the system (cos fkx)kLi is a Riesz basis of Ti 
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(see Lemma [A. 1|) and the numbers := — cos/fc = (— l) fc sin define an £ 2 -sequence. 
Therefore there exists a unique function g 6 7i with Fourier coefficients in the Riesz 
basis (cos fkx)^ =1 . Now the function 

(3.6) /(A):=cosA + / g{x) cos Ax dx 

Jo 

is an even entire function of order 1 and has zeros at the points ±fk- It follows that / 
and / differ by a scalar factor, which as before is shown to equal 1. 

Due to the asymptotics of the zeros fk only finitely many of them can repeat. If, 
e.g., /„ = f n+ \ = ■ ■ ■ = f n+p -i is a zero of / of multiplicity p = p(n), then we include 
to the above system the functions cos f n x, x sin f n x, ■ ■ ■ , x p ~ l sin(7rp/2 — f n x). After 
this modification has been done for every multiple root, we again get a Riesz basis 
of H and can find a function g 6 7i with Fourier coefficients c n := — cos/ n ,c„ + i = 
— sin /„,•■■ , c n+p („) := — sin(7rp(n)/2 — f n x). Then ±f n becomes a zero of the function 

/ in ()3.6|) of multiplicity p(n), and the proof is completed as above. □ 
Lemma 3.3. The following equalities hold: 

$ 1 (A) = n^ 2 (^-V2)- 2 (A 2 fe -A 2 ), 
3 ' 7 $ 2 (A) = n^ 2 (^-V2)- 2 (^-A 2 ). 

Proof. We recall ^2] that the function u\ has a representation via a transformation 
operator connected with the point x — 1 of the form 

-i 



Ui(x, A) = cos A(l — x) + / k(x, t) cos A(l — t) dt, 
which implies that 



$i(A) = cosA+ / fci(t) cos Xtdt, 
Jo 



where k\{t) := k(0, 1 — t) is a function from Ti. The claim now follows from Lemma[ 
The representation for $2 is derived analogously, and the lemma is proved. □ 

Now, given two spectra (A 2 ) and ) of two Sturm-Liouville operators T a ^ hl and T a ^ h2 
respectively, with unknown a e H and hi,h 2 G R, we proceed as follows. First, we 
identify hi — h 2 as lim n ^ 00 (A 2 — fil)/2 (see Corollary 12. 4|) . then construct the functions 
$1 and $2 of ()3.7|) and determine the norming coefficients (ct„) via 1)3.1)) . The spectral 
data {(A 2 ), (a n )} determine a and /ii up to an additive constant by means of the 
algorithm of [llj. This gives the required function a and two numbers hi and h 2 up to 
an additive constant, and the reconstruction is complete. 

The second part of the inverse spectral problem is to identify those pairs of sequences 
(A 2 ) and (/i 2 ) that are spectra of Sturm-Liouville operators T ahx and T ah2 for some 
real- valued a G H and some real hi,h 2 , hi 7^ h 2 . We established in Section |21 the 
necessary conditions on the two spectra given by Theorem 12.11 and Lemmata 12.21 and 
12.3) it turns out that these conditions are sufficient as well. Namely, the following 
statement holds true. 

Theorem 3.4. Suppose that sequences (A 2 ) and (/i 2 ) of positive pairwise distinct num- 
bers satisfy the following assumptions: 

(1) the sequences (A 2 ) and (/i 2 ) interlace; 
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(2) A n = n(n — 1/2) + X n and \i n = n(n — 1/2) + fi n with some £ 2 -sequences (A n ) 
and (/I n ); 

(3) there exist a real number h and an t2-sequence (z/„) such that X n — fi n = + 

Then there exist a function a G H and two real numbers hi and h 2 such that (A n ) 
and (fi n ) are the spectra of the Sturm- Liouville operators T a hl and T ajh2 respectively. 
Moreover, the operators T a ^ x and T a ^ 2 are recovered uniquely. 

In order to prove the theorem we have to show first that the numbers a n constructed 
for the two sequences through formula (|3.1j) . with the functions $i,$2 of (|3.7j) and 
h\ — h 2 := h with h of item (3), are positive and obey the asymptotics a n = 1 + a n 
for some ^-sequence (5 n ). Then the algorithm of [TT] uses (A^) and (a n ) to determine 
(unique up to an additive constant) function a G Ti and hi G M such that {A^} is the 
spectrum of the Sturm-Liouville operator T (Ti / ll and a n are the corresponding norming 
constants. The second and final step is to verify that {/i^} is the spectrum of the 
Sturm-Liouville operator T a ^ 2 with h 2 '■— hi — h, where h is the number of item (3). 
These two steps are performed in the following two lemmata. 



Lemma 3.5. Assume (l)-(3) of Theorem \3-4\ Then the numbers a n constructed 
through relation (|3.1j) with $i,$2 of (|3.7|) and hi — h 2 := h with h of item (3) are 
all positive and obey the asymptotics a n = 1 + a n , where (a n ) is an i 2 -sequence. 

Proof. By Lemma l3~2l there exist 7i-functions f\, f 2 such that the functions $j, j = 1, 2, 
of ()3.7|) admit the representation 



$j ( A) = cos A + / fj (t) cos At dt. 
Jo 



It follows that 



(A n ) = - sin A„ - / tf\ (t) sin X n t dt 
Jo 



= (-l) n cosA n - / tfi{t)smX n tdt = {-l) n + X n , 
Jo 

where (X n )^Li G £ 2 , and by similar arguments 

$ 2 (A n ) = (A n - /i n )$ 2 (/i„) + 0(\X n - IJ n \ 2 ) 

= (A n -/i n )[(-ir + /y 

for some ^2-sequence (fin)- Next, assumptions (2) and (3) easily imply the relation 

h 1 - 

= 1 + v n , 



A n (A n — A 4 

for some ^-sequence (0 n ). Therefore the numbers 

_ h 

Otri. 



X n $ 2 (A r 



obey the required asymptotics. 

Finally, the interlacing property of the two sequences implies that all a n are of the 
same sign, and thus are all positive in view of the asymptotics established. The proof 
is complete. □ 
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With the above-stated properties of the sequences (A^) and (a n ) we can employ the 
result of [TT] that guarantees existence of a unique Sturm-Liouville operator T a>hl with 
a real- valued function a G Ti and a real number hi such that {A 2 } coincides with 
the spectrum of T a>hl and a n are the corresponding norming constants. Now we put 
hi := h\ — h with h of assumption (3) and expect {/i 2 } to be the spectrum of T a ^ 2 - 

Lemma 3.6. The spectrum of the Sturm-Liouville operator T a ^ 2 with the above a Eli. 
and h 2 G M coincides with (/x 2 ). 

Proof. Denote by v i(x, A) a solution of the equation l a {y) = X 2 v satisfying the boundary 
conditions t>(l) = l,t>M(l) = —hi. Then Vi(x, A n ) is an eigenfunction of the operator 
T a)hl corresponding to the eigenvalue A 2 , and by construction 

2/ |wi (x, A n )| dx = a n := — 

JO ^n < ^2{^n) 

with the functions $i, $ 2 of (|3~7|) . 

Denote by (i^) the spectrum of T a ^ 2 an d put 

oo 

$ 2 (A) :=n^ 2 (^-V2)- 2 (^-A 2 ). 
fc=i 

By Lemma f3. II we have 

$UA„) 
K $ 2 (A n ) 

whence $ 2 (A n ) = $2(A n ) for all neN. Lemma f3.2l implies that 

$ 2 (A) - $ 2 (A) = / f(t)cos\tdt 
Jo 

for some f E H. Equality $ 2 (A n ) = $2(A n ) means that the function / is orthogonal to 
cos X n t, n6N. Since the system {cos \ n t} forms a Riesz basis of TC, we get / = and 
\P 2 = \l/ 2 . Thus = fi n , and the lemma is proved. □ 



4. Reconstruction by Dirichlet and Dirichlet-Neumann spectra 

The analysis of the previous two sections does not cover the case where one of hi, hi 
is infinite. In this case the other number may be taken without loss of generality 
(recall that T a)h = T a+h)0 ), i.e., the boundary conditions under considerations become 
Dirichlet and Dirichlet-Neumann ones. 

Suppose therefore that a G TC is real valued and that (A 2 ) and (/i 2 ) are spectra of 
the operators T a>0 and T a>00 respectively; without loss of generality we assume that A n 
and n n are positive and strictly increase with n. The reconstruction procedure remains 
the same as before; namely, we use the two spectra to determine a sequence of norming 
constants (a n ) and then recover a by the spectral data {(A 2 ), (a n )}. 

Denote by w_(-,A) and «+(•, A) solutions of the equation l a u = \u satisfying the 
initial conditions w_(0, A) = 0,m^(0, A) = A and w + (l,A) = l,tt+(l,A) = respec- 
tively. Then « + (-,A n ) is an eigenfunction of the operator T aj0 corresponding to the 
eigenvalue A 2 and 

a n := 2 / \u+(x, A n )| 2 dx 
Jo 
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is the corresponding norming constant. We also put 1 J r i(A) := u+(0, A) and ^(A) := 
u_(l, A); then zeros of the functions and ^2 are numbers ±A n and ±/i n respectively. 
As earlier, the function ^ is uniquely determined by its zeros through formula (J3.7j) . 
Observe that the Dirichlet eigenvalues /z 2 have asymptotics different from that of A 2 , 
so that ^2 requires slight modification of formula f)3.7j) . We first investigate the asymp- 
totics of fl n . 

Theorem 4.1 (^21 EI)- Suppose that o e 7i and {^}, n GN, is the spectrum of the 
operator T a}00 . Then the numbers fi n satisfy the relation 

fx n = 7m + Jl n , 

in which the sequence (jj, n ) belongs to £2- 

Proof. The solution U- can be represented by means of the transformation operator ^2] 
as U-(x, A) = sin Aa; + f Q k£(x, t) sin Xt dt, where k+ is the kernel of the transformation 
operator. Therefore the numbers ±/i„ are zeros of the function 

(4.1) * 2 (A) := sinA + / fc+(l, t) sin At dt, 



which is entire of order 1. Since the function k£(l, •) belongs to Ji the required 
asymptotics of \i n is derived in a standard way (cf. [TH] and Section |2J . □ 

The function ^2 is determined by its its zeros in the following way. 

Lemma 4.2. The following equality holds: 

(4.2) * 2 (A) = An^~V fe -A 2 ). 

The proof is completely analogous to that of Lemma 13.21 and is left to the reader. 
Now we show how the norming constants a n are expressed via ^1 and ^> 2 - 

Lemma 4.3. The norming constants satisfy the following equality: 

Proof (cf. [5\). The Green function G(x,y, A 2 ) of the operator T^o (i. e., the kernel of 
the resovlent (T^o — A 2 ) -1 ) equals 

nl X 2n 2 u + (x,X n )u + (y,X n ) 
G(x,y,X) = }_^——- 

On the other hand, we have 



G(x,y,X 



2\ 



U-(x, X)u + (y, A), < x < y < 1, 
W(X) '\u+(x,X)u4y,X), Q<y<x<l, 

where W(X) := u+(x, X)u m (x,X) — U-(x, A)w+ (x, A) is the Wronskian of u + and 
The value of VT(A) is independent of x G [0, 1]; in particular, taking x = we find that 
W(X) = Xu+(0,X) = A*i(A). 

Now we take x = y = 1 in the above expressions and find that 

y> 2 1 = tt 2 (A) 

^A 2 -A 2 « n - A^(A)' 

Comparing the residues at the poles A = A n , we derive formula (|4.3|) . and the lemma 
is proved. □ 
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Now the reconstruction procedure is completed as follows: we determine the sequence 
of norming constants (a„) as explained above and find a unique operator T a0 with 
spectrum (A^) and norming constants (a n ), see [UJ. This gives the function a G Ti 
and thus the operator T a oo . 

Now we would like to give an explicit description of the set of all possible spectra 
of the operators T CT) o and T a>00 when a real- valued function a runs through Ti, i. e., 
to give the necessary and sufficient conditions on two sequences (A^) and (n n ) to be 
the spectra of T^o and T CT)00 with a real-valued a G Ti. Necessary conditions are that 
the eigenvalues X n and n n should interlace and obey the asymptotics described in 
Theorems 12. II and 14. 1\ we shall show that these conditions are in fact sufficient as well. 

Theorem 4.4. Suppose that sequences (X n ) and (fi n ) of positive pairwise distinct num- 
bers satisfy the following assumptions: 

(1) the sequences (X n ) and (fi n ) interlace: < Af < //f < A| < . . . ; 

(2) X n = ir(n — 1/2) + X n and fi n = 7in + Jl n with some l 2 -sequences (A n ) and (fin)- 

Then there exist a unique function o G Ti such that (X n ) and (fi n ) are the spectra of 
the Sturm- Liouville operators T CT) o and T CTi00 respectively. 

The main ingredients of the proof of this theorem are the same as for Theorem 13.41 
first, we construct functions and ^/ 2 by their zeros and define the sequence (a n ) 
as explained in Lemma I4.3| then we prove that a n are all positive and have the re- 
quired asymptotics ot n = 1 + a n for some ^-sequence (d n ). Using now the reconstruc- 
tion algorithm of [TI], we find a unique function o G Ti such that the correspond- 
ing Sturm-Liouville operator T ct0 with potential q = a' G Ti possesses the spectral 
data {(A^), (a n )}. Finally, we prove that (fi n ) is the spectrum of the operator T G)00 
with a just found, and the reconstruction procedure is complete. 

5. The case of the Neumann boundary condition at x = 

The analysis of the previous sections can easily be modified to cover the case H = 0, 
i.e., the Neumann boundary condition iiW(0) = 0. As before, we use the two spectra 
to determine the sequence of norming constants and then apply the reconstruction 
procedure of [TTJ to find the corresponding Sturm-Liouville operators. Also the neces- 
sary and sufficient conditions on the two spectra can be established. We formulate the 
corresponding results in the following two theorems. 

Theorem 5.1. Suppose that sequences (X n ) and (fi n ) of positive pairwise distinct num- 
bers satisfy the following assumptions: 

(1) the sequences (X n ) and (fi n ) interlace; 

(2) A„ = n(n — 1) + A„ and fi n = n(n — 1) + Jl n with some ^-sequences (A n ) and 
(fin); 

(3) there exist a real number h and an £ 2 sequence {y n ) such that \ n — [i n = ■^ + i n L - 

Then there exist a unique function a ETi and real constants hi, h 2 such that (\ n ) and 
(fi n ) are the spectra of the Sturm-Liouville operators T a o thl and T a o h2 respectively. 

Conversely, the spectra (X n ) and (/x^) of Sturm-Liouville operators T^o,^ and T at0 ^ 2 
with a <ETi and hi,h 2 6l satisfy assumptions (l)-(3). 

For the case where one of hi,h 2 is infinite (say, h 2 = 00) the asymptotics of the 
corresponding spectrum is different; also assumption (3) becomes meaningless and 
should be omitted. 
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Theorem 5.2. Suppose that sequences (A^) and (/i^) of positive pairwise distinct num- 
bers satisfy the following assumptions: 

(1) the sequences (A^) and (/i^) interlace; 

(2) A„ = ir(n — 1) + A ra and /i n = n(n — 1/2) + /I n with some ^sequences (A n ) and 

Ov)- 

ITien t/iere exist a unique function a G and a reaZ constant hi such that (A^) and 
(/i^) are tne spectra of the Sturm- Liouville operators T^o^ and T^oo respectively. 

Conversely, the spectra (A n ) and (/i^) o/ Sturm- Liouville operators T^o,^ and T^oo 
u>ztn a Eli. and h\ G R. satisfy assumptions (1) and (2). 

Appendix A. Riesz bases 

In this appendix we gather some well known facts about Riesz bases of sines and 
cosines (see, e.g., [HI IE] and references therein for a detailed exposition of this topic). 

Recall that a sequence (e n )J° in a Hilbert space Ti is a Riesz basis if and only if 
any element e 6 H has a unique expansion e = Y^=i c n e n with (c n ) G £2- If (e n ) 
is a Riesz basis, then in the above expansion the Fourier coefficients c n are given by 
c n = (e, e' n ), where (e^)f 3 is a system biorthogonal to (e n ), i.e., a system which satisfies 
the equalities (efe,e^) = 5k n for all k,n G N. Moreover, the biorthogonal system (e^) is 
a Riesz basis of Ti as long as (e n ) is, in which case for any e G 7i also the expansion 
e = XX e ' e ™) e n takes place. In particular, if (e n ) is a Riesz basis, then for any e e H 
the sequence (c' n ) with c n := (e, e n ) belongs to £ 2 - 

Proposition A.l ([E])- Suppose that Hk — > as A; — > 00 and taat tae sequence nh + fik 
is strictly increasing. Then each of the following systems forms a Riesz basis o/L 2 (0, 1): 

(a) {sin(7rfcc + jj, k x)}™ =l ; 

(b) {sin(7r[A;-l/2]x + /x fc x)}^ =1 ; 

(c) {cos(irkx + fi k x)}^L ; 

(d) {cos(vr[A; + l/2]x + /i fc x)}^ =0 . 

Lemma A. 2. Suppose that (A n ) and (/i n ) are two sequences of real numbers such that 
lim [A n - 7r(n - 1/2)1 = lim \/i n - ir(n - 1/2)1 = 0. 

n — >oo ra— >oo 

and assume that f E TC. Then there exists an t2-sequence [y n ) such that 
1 

/(t)[cosA n t - COS /i„t] dt = (A n - \ln)V n + 0(|A„ - /i n | 3 ) 
as n — > 00. 

Proof. Using the relation 

cosA n t - cos fi n t = -2sin[(A n - //„)t/2] sin[(A n + fi n )t/2], 

we find that 

/ /(t)[cosA n t - cos/i n t] dt = -2 / f(t) sin[(A„ - ix n )t/2] sin[(A n + //„)t/2] dt 
Jo Jo 

= -(A n -/i n ) / t/(t)sin[(A n + /^t/2]dt + 0(|A n -/^| 3 ). 
Jo 

We put v n := — J Q tf(t) sin[(A n +// n )t/2] dt and observe that the sequence (^„) belongs 
to £2 since the system {sin[(A n + /i n )t/2]} is Riesz basic in Ti (at least for all n large 
enough) and the function tf(t) belongs to Ti. □ 
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